Constructing a Large Variety of Dirac-Cone Materials in the Bii-j^Sbj^ Thin Film System 
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We theoretically predict that a large variety of Dirac-cone materials can be constructed in Bii_;(;Sbx thin films, 
and we here show how to construct single-, bi- and tri- Dirac-cone materials with various amounts of wave vector 
anisotropy. These different types of Dirac cones can be of special interest to electronic devices design, quantum 
electrodynamics and other fields. 
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Dirac cone materials have recently attracted considerable 
attention. In an electronic band structure, if the dispersion re- 
lation E{k) can be described by a linear function as £^ = v • k, 
where v is the velocity, k is the lattice momentum, and h= 1, 
the point where £^ ^ is called a Dirac point. A Dirac cone 
is a two-dimensional (2D) Dirac point. Dirac cone materials 
are interesting in electronic device design, quantum electro- 
dynamics and desktop relativistic particle experiments etc. A 
single-, bi- or tri-Dirac cone system has one, two or three dif- 
ferent Dirac cones degenerate in E(k) in the first Brillouin 
zone. Graphene has two degenerate isotropic Dirac cones at 
points K and K' in its first Brillouin zone, which is there- 
fore considered as a bi-Dirac-cone system. Many novel phe- 
nomena are observed in this system [1], such as the room 
temperature anomalous integer quantum Hall effect the 
Klein paradox which means that fermions around a Dirac 
cone can transmit through a classically forbidden region with 
a probability of 1. Dirac fermions can be immune to local- 
ization effects and can propagate without scattering over large 
distances on the order of micrometers ||4i] . 

In this Letter, we show how to obtain single-, bi- and tri- 
Dirac-cone Bii__;(:Sbx thin films, and how to construct Dirac 
cones with different anisotropics. We also point out the possi- 
bility of constructing semi-Dirac cones in Bii_;cSbx thin films. 

Bii_;cSb;c has many special properties that are interesting 
from the point of view of anisotropic Dirac cones. We re- 
call that bulk Bii_;cSbx is a crystalline alloy with a rhombo- 
hedral structure, which displays remarkable anisotropy. The 
first Brillouin zone of bulk Bii_;(:Sbx has one T point and three 
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FIG. 1: The 3-fold degenerate L-point electron pockets and the T- 
point hole pocket in the first Brillouin zone of bulk Bii-j^Sbx. 



degenerate L points, L^^\ L^^^ and L^^ , as illustrated in Fig. 1 . 
The bottom of the conduction band is located at the L points, 
while the top of the valence band can be located either at the 
T point or at the L points, depending on the Sb composition 
X when < x < 0.10. In bulk Bii__;cSbx, the band structure 
varies as a function of Sb compostion x, temperature T, pres- 
sure P and stress t f?]. The conduction band is very close 
to the valence band at the L points, so that these bands are 
non-parabolically dispersed as [6] 



£(k) = ±((vk)2 + £:; 



(1) 



due to their strong interband coupling. When the L-point band 
gap Eg is small, the dispersion relation E{k) becomes linear 
and Dirac points are formed as E{k) ^ ±v • k. The L-point 
band gap Eg can approach under some conditions, e.g. when 
P = 1 atm. Eg ^Ositx^ 0.04 and T < 11 K ov Sit x ^ 
0.02 and L < 300 K [8]. For simplification, this Letter will 
focus on the low temperature range (T < 11 K) where the band 
structure of Bii-j^Sbj^ does not change much with temperature. 

For Bii_;tSbx thin films, the 2D band structure varies also as 
a function of film thickness and film growth orientation, which 
provides considerable flexibility compared to bulk Bii_;cSbx. 
Furthermore, the quantum confinement effect in the thin film 
system is potentially interesting, where its anisotropic proper- 
ties imply potential application possibilities. 

The energy spectrum near an L-point Dirac cone in a 
Bii_xSbx thin film is calculated based on the iterative-two- 
dimensional-two-band model described below. Here the gen- 
eral two band model for two strongly coupled bands obeys the 
relation |0] 



p-a-p = £(k)(i + :^), 



(2) 



where p is the carrier momentum vector and 06 is the inverse- 
mass tensor. The two coupled key parameters a and Eg are 
calculated in an iterative way in our model as 
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FIG. 2: An illustration of (a) single-, (b) bi- and (c) tri-Dirac-cone 
Bii_;cSbjc thin films grown along the (a) bisectrix, (b) binary and (c) 
trigonal axes, respectively. For the cross-sectional view of each cone, 
k is chosen such that V\^E(k) has its minimum along that direction 
of k. The illustration is based on the example of Bii_;(-Sbx thin films 
with Iz = 100 nm, x = 0.04, P = 1 atm and T < 77 K, under which the 
L points of bulk Bii_;^-Sbx have a zero-gap. The scenario is similar 
for other conditions. In (a), a single-Dirac-cone is formed at the L^^"^ 
point, while the L*^^^- and L*^^^- point band gaps are opened up. In 
(b), two degenerate quasi-Dirac cones are formed at the L*^^) and L*^^) 
points, while the L*^^^ -point band gap is much larger, which leads to 
a bi-quasi-Dirac-cone material. The band gap at the L^^'^ and L^^'^ 
points can be less than 1 meV if a sample of 1^ = 200 nm is chosen, 
which leads to exact Dirac cones. In (c), the L^^)- and L^^)- 
point band gaps are all the same, and the three quasi-Dirac cones are 
degenerate in energy. 
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where mo is the free electron mass, I is the identity matrix, 
Iz is the film thickness and n denotes the nth step in the it- 
eration. The procedure is repeated until a'^l and E^^^ be- 
come self-consistent, and then we get accurate solutions for 



af'^^{Bii_^Sbj,) = aW and E^'^'^iBiu 

Bii_;tSbx. Because of the af 
for Sb composition < x < 0.10, 



g ^ for thin 

film Bii_;tSbx. Because of the approximations that are valid 
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FIG. 3: The anisotropy of the L*^^) -point Dirac cone vs. film growth 
orientation, (a) The anisotropy coefficient 7 vs. film growth orienta- 
tion. The value of /for a specific film growth orientation is shown by 
the radius and color. 7 can be as large as ~14 for films grown along 
the trigonal axis, and as small as ~2 for films grown along the bisec- 
trix axis. V is shown for the L*^^) -point Dirac fermions vs. transport 
direction for Bii-^^Sbx thin films grown along the (b) trigonal, (c) 
bisectrix and (d) binary axes, (b)-(c) are drawn based on an example 
sample with 1^ = 300 nm and x = 0.04. 



Eqs. ([3]) and dH) can be further simplified, which converge to 
the analytical solution as 



a' 



bulk 



E^/'\Bi) (5) 



and 



7t^ a{^""{Bh.,Sb,) 



Ef"'"iBh.M = E'/'iBh.M - 

(6) 

The dispersion relation E{k) can then be the solved by the 
methods used by Ref . from 
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where OLij = a/sa^s/ass — atj for ij = 1 and 2, and a = 
af'^^{Bii_^Sb^). The Hamiltonian for Bi and Bii_;,Sb;, based 
on k • p theory in Eq. (O is equivalent to a Dirac Hamilto- 
nian with a scaled canonical conjugate momentum 11 1111 . Thus, 
Eq. ©is also a good approximation to describe the Dirac 
cones. The band parameters we use in the present calcula- 
tions are values that were measured by cyclotron resonance 
experiments fl5]. 

According to Eqs. ([B and (|7]), when £'g ^ at an L point, 
the electronic dispersion relation becomes a perfect Dirac 
cone, where the energy E is exactly proportional to the lattice 
momentum k measured from that L point. When Eg becomes 
large enough lO . the linearity of the dispersion relation be- 
comes an approximation, and the Dirac cone becomes a quasi- 
Dirac cone. If an ^ 0.22 with a finite Eg, so that E oc kx and 
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FIG. 4: Scheme for the L*^^) -point band gap vs. film growth orienta- 
tion, film thickness and Sb composition. (a)Illustration of the L*^^)- 
point band gap vs. film growth orientation and film thickness. The ra- 
dius, direction and color represent the film thickness, film growth ori- 
entation and L*^^) -point band gap, respectively. The illustration takes 
X = 0.04 as an example. For other Sb compositions (0 <x < 0.10), 
the film thickness and film growth orientation dependence for the 
L^^) -point band gap should be similar, which is here illustrated for 
thin films grown along the (b) trigonal, (c) bisectrix and (d) binary 
axes. 
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ken. The value of the film-direction-inverse-mass-component 
ay^^{Bii_xSbx) is much smaller for the L^^^ point than the 
corresponding values for the L^^^ and L^^^ points. The L^^^- 
point gap E^^ is negligibly small due to the small value of 
a^2^{Bii_xSbx), where a Dirac cone is formed, as shown in 
Fig. 2a. However, the L^^^- and L^^^- point band gaps Eg^'' and 

(3) 

Eg are much larger, which implies that a single-Dirac-cone 
at the L^^) point is constructed. Here we are taking advantage 
of both the extreme anisotropy of Bii-j^Sb^ and the quantum 
confinement effect of thin films. The quantum confinement ef- 
fects for the L^^) -point carriers differs remarkably from those 
for the L^^)- and L^^^ -point carriers due to the anisotropy of 
the L-point pockets. Figure 2b shows that a Bii-^Sbjc thin 
film grown along the binary axis can be a bi-Dirac-cone mate- 
rial, where the L^^^ -point band gap Eg^'' is much larger than the 
L^^)- and L^^^- point band gaps Eg'^'' and Eg^\ Thus, £'g^^ and 

(3) 

Eg ^ remain small enough [13] to make two degenerate Dirac 
cones (quasi-Dirac cones) at the L^^^ and the L^^^ points. In 
Fig. 2c, the film is grown along the trigonal axis, so that 
the 3 -fold symmetry for the three L points is retained. The 
three Dirac cones (quasi-Dirac cones) at the L^^\ L^^^ and L^^^ 
points are degenerate in energy, which makes this film a tri- 
Dirac-cone material. By definition, an exact Dirac cone has 
Eg = 0. However, Eg = Q Dirac cones are seldom achieved 
experimentally, so it is practical to consider £'g < /^gT as a 
criterion for an exact Dirac cone. In the temperature range 
below 77 K that we are considering in this paper, the thermal 
smearing of ksT corresponds to ~ 7 meV. For the criterion of 
a quasi-Dirac cone, we can use ksT <Eg< Eg{Bi)^^^^, where 
EgiBif""^^ 14 meV. Thus, we consider the three Dirac cones 
in Fig. 2c, as quasi-Dirac cones, which are plotted for the 
case of = 100 nm and EgC::i\Q meV. If exact Dirac cones 
are needed, a larger film thickness can be chosen, e.g. = 200 
nm, which satisfies Eg <kBT . 

We now show how to construct anisotropic Dirac cones 
with different shapes for the wave vector as a function of cone 
angle. To characterize the anisotropy of a Dirac cone, we de- 
fine an anisotropy coefficient 



FIG. 5: Example of a semi-Dirac cone in the Bii-j^Sbx thin film 
system (x = 0.10 and = 100 nm). It can be seen that around the 
L^^^ point, the fermions are relativistic (linearly dispersed) along the 
Vmax direction, and classical (parabolically dispersed) along the Vmm 
direction. 



E oc ky, we call it a semi-Dirac cone. In a semi-Dirac cone, 
the fermions are relativistically dispersed in one direction (kj), 
and classically dispersed in another direction (ky). 

We propose that single-, bi- and tri-Dirac-cone materials 
can be constructed from Bii-^^Sb^ thin films, by proper syn- 
thesis conditions to control the relative symmetries of the 3 L 
points. Bii_;cSbx thin films grown along the bisectrix axis can 
be single-Dirac-cone materials, as illustrated in Fig. 2a, where 
the 3-fold degeneracy of the L^^\ L^^^ and L^^^ points is bro- 
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where Vmax and Vmin are the maximum and minimum in-film 
carrier group velocities for a Dirac cone that is defined as 
v(k) = Vk£'(k). For a perfect Dirac cone, v is a function of 
the direction of the lattice momentum k measured from that 
L point only and is independent of the magnitude of k. For 
an imperfect Dirac cone or a quasi-Dirac cone, this magni- 
tude invariance is exact only when k is large, and becomes an 
approximation around the apex when k is small [13]. 

Fig. 3 gives us an important guide on how to construct 
anisotropic L^^^ -point Dirac cones. In Fig. 3a, the anisotropy 
coefficient y for the L^^^ -point Dirac cone as a function of film 
growth orientation is shown. For a film grown along the bi- 
sectrix axis, 7 has its minimum value Ymin = "^2, where the 
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carrier velocity v(k) for the L^^^ -point Dirac cone varies only 
by a small amount with the direction of k, as shown in Fig. 3b. 
For a film grown along the binary axis, 7= ~10, where v(k) 
varies more with the direction of k as shown in Fig. 3c, com- 
pared to Fig. 3b. For a film grown along the trigonal axis, 7 
has its maximum of Jmax = ~14, where v varies significantly 
with the direction of k, as shown in Fig. 3d. 

Researchers have tried to realize semi-Dirac cones in oxide 
layers lfl4ll . where the fermions are relativistic in one direction 
and classical in its orthogonal direction. In the present work, 
we have found that it is possible to construct semi-Dirac cones 
in the Bii_;cSbx thin film system. According to Eqs. ©, © 
and (|7]), for an in-film direction k, where k is a unit directional 
vector of k in the in-film lattice momentum space, whether 
the dispersion relation is linear or parabolic depends on the 
L-point band gap Eg, and the a projection along that direction 
of k, defined by OJj^ = k* • a • k, where d is given by Eq. (|7]). 
When Eg is small and is large lll3ll . the energy becomes 
linearly dispersed along k; when Eg is large and is small, 
the energy becomes parabolically dispersed along k. To con- 
struct a semi-Dirac cone, we needs to find a proper L-point 
band gap Eg and anisotropy 7, such that Lg/Omax is small 
and Lg/Omin is large. In this case, the electronic energy is 
linearly dispersed along the Omax direction and parabolically 
dispersed along the Omin direction. Here o^max and Omin are 
the maximum and minimum values of a^, which correspond 
to the principal axes of the 2D tensor d. The L-point band 
gap varies as a function of the film thickness Z^, the growth 
orientation and the Sb composition x, as shown by the calcu- 
lated results given in Fig. 4. To construct a semi-Dirac cone, 
we need to find a growth direction that ensures a significant 
anisotropy, and a large enough value of Eg which ensures that 
L(k) becomes parabolically dispersed along the Omin direc- 
tion. However, the Eg should not be too large, because of the 
necessity that the linear dispersion relation along the Omax di- 
rection is maintained. These requirements can all be achieved 
by choosing the proper Sb composition x, film thickness 1^ and 
growth orientation as shown in Fig. 4. From Figs. 2 and 3, we 
know that the L^^^ -point Dirac cone has a maximum k- vector 
anisotropy, when the growth orientation is near the trigonal 
axis. We also see that the thin film becomes a single-Dirac- 
cone material when the growth direction is near the bisectrix 
axis. Thus, a good strategy to construct a semi-Dirac cone is 
to choose a growth orientation between the trigonal and the 
bisectrix axis, in the trigonal-bisectrix plane. An example of a 
semi-Dirac cone is shown in Fig. 5, where the example sam- 
ple is grown along a direction that is 40° from the trigonal 



axis, 50° from the bisectrix axis, and perpendicular to the bi- 
nary axis. Thus, a large Sb composition (e.g. x ^ 0.10) and a 
small film thickness (e.g. ^ 100 nm) is preferred to make 
the Eg large, and x = 0. 10 and = 100 nm are chosen for this 
example sample. 

In conclusion, we have proposed the growth of Bii_;cSbx 
thin films, which for selected concentrations of Sb and differ- 
ent directions to the film normal allows different Dirac-cone 
materials to be constructed. We have shown how to construct 
single-, bi- and tri-Dirac-cone materials, as shown in Fig. 2, 
as well as quasi- and semi-Dirac-cone materials, as shown in 
Fig. 2c and Fig. 5, respectively. 
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